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ABSTRACT 


Laser-Doppler velocimeter measurements were made in the flow field of a stalled 
cascade of controlled-diffusion stator blades. Tests were conducted at 10 degrees of 
incidence above the design inlet-flow angle in order to verify previous measurements. A 
unique method of data presentation was offered to characterize the unsteady positive and 
negative velocities about their mean. Laser-sheet flow visualization was performed and 
showed the unsteadiness of the stalling phenomena within the blade row passages. 
Additionally, a leading-edge separation bubble was observed and reversed flow was 
measured within the bubble using the laser-Doppler anemometer. Power-spectral density 
and auto-correlation analysis results are presented. In parallel with the experimental 
measurements, a computational fluid-dynamics study was initiated in an attempt ultimately 
to predict stall. Viscous flow computations were completed at on-design and near-stall 
conditions. Pressure distributions, separation bubble re-attachment locations, and loss 
predictions were compared to previous experimental studies. Generally good agreement 


was obtained for on-design conditions with degradation in prediction nearing stall. 
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I INTRODUCTION 


Compressor stall and off-design behavior impose performance constraints on gas 
turbine engines. The open literature currently lacks adequate documentation of 
quantitative measurements of turbomachinery flow fields at and beyond stall. This 
knowledge is essential for predicting the onset of stall with computational methods and 
optimizing blade geometries for off-design range with minimal degradation in design-point 
performance. Towards optimum design-point performance, controlled-diffusion (CD) 
blading was developed. The purpose of these blades was to avoid boundary layer 
separation and, in the case of supercritical (SC) blading, to be shock-free in the transonic 
range. 

Extensive laser-Doppler velocimeter (LDV) investigations [Refs. 1 - 3] of flow 
through a controlled-diffusion compressor cascade have been conducted in the region 
below stall at the Naval Postgraduate School (NPS). Recently, the reported cascade was 
stalled with an inlet-flow angle of 50 degrees and an extensive mapping of the flow field 
was performed using a two-component LDV system by Ganaim [Ref. 4]. The motivation 
for the present study was to partially map the flow field near stall and to establish the 
repeatability of the results reported in Reference 4. Between Ganaim’s work and the 
present, the tunnel was reconfigured twice to pre-stall flow angles. This was done to 
confirm the repeatability of the flow field and to expand the available database for 
unsteady separated flows. Also, since the previous data were presented on a mean-flow 
basis, it was desired to more properly characterize the mean flow into its positive and 
negative components. This required a special method of post-processing the velocity- 
histogram data, as described within this report. 

It was also desired to characterize the unsteady flow using power-spectrum and 
auto-correlation analyses. It was hypothesized that there could exist a fundamental 
frequency for the stalled cascade, that it could be determined from the measurements made 
in the present study and should match that of the previous data [Ref. 4]. 

In parallel with the LDV study, a computational fluid dynamics (CFD) study was 


initiated in an attempt ultimately to predict stall. Previous CFD studies on the present 


cascade include those of Elazar [Ref. 1], Stow, et al [Ref. 5], Hobson [Ref. 6], and Kang, 
et al [Ref. 7]. 

Before attempting to calculate the unsteady, stalled, flow through a cascade (as 
was measured in the present study), it was necessary first to evaluate the code’s predictive 
ability at incidence angles off-design and near the stall point. A quasi-three-dimensional 
computer code [Refs. 8-9] for the solution of the thin-layer Navier-Stokes equations was 
used to predict the cascade performance from design (40 degrees inlet-flow angle) to near 
stall (48 degrees inlet-flow angle). The code was used to successfully compute the blade 
surface pressure distribution, the leading edge separation bubble and total pressure loss 


coefficient, for inlet flow angles from 40 to 48 degrees. 


Il. EXPERIMENTAL SETUP 


A. LOW-SPEED CASCADE WIND TUNNEL 


The subsonic cascade wind tunnel contained 20 controlled-diffusion blades. Each 
blade was 10 inches in span, 5.01 inches in chord and separated by 3 inches in the 
pitchwise direction. The cascade solidity (blade chord-to-pitch ratio) was 1.67 and 
Reynolds number based on chord was approximately 700,000. A detailed description of 
the test facility, test section and CD blading was fully documented in Reference 1. Figure 
1 shows a schematic of the cascade. The tunnel was configured for an inlet-flow angle 
(B:) of 50 degrees, and an exit-flow angle ($2) such that the downstream end-wall static 


pressure distribution was uniform. 


B. INSTRUMENTATION 


Blade surface pressure measurements were recorded using a 48-channel Scanivalve 
system controlled by an HP-9000 computer. Figure 1 shows the “fully-instrumented 
blade” location within the test section. Reference 1 gives the location of the blade surface 
pressure points. A full description of the data acquisition system is given in Reference 10. 

A two-component, four-beam TSI Model 9100-7 LDV system was used for all 
flow field measurements. The LDV setup, including laser type, optics, atomizer, seeding, 
traverse mechanism and data acquisition were thoroughly documented in Reference 4. 

The flow was visualized by injecting fog particles of 2 um in size through a wand 
inserted through the end wall upstream of the blading. A 300 mW air-cooled argon-ion 
laser was used as the light source for laser-sheet flow visualization. The laser beam was 
transmitted to a lens system via a fiber-optic cable. The lens system ended with a 
cylindrical-plano lens, which produced the laser sheet. The laser sheet was introduced into 
the cascade tunnel, ahead of the blades, through the sidewall. A schematic of the system, 
showing the location of the introduction of the fog and the position of the video camera 


(located outside the window), are shown in Figure 2. 
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Figure 1. Schematic of the Cascade Wind Tunnel 
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C: EXPERIMENTAL PROCEDURE 


1. Surface Pressure Measurement 


The tunnel plenum chamber pressure was set at 12 inches of water (gauge) and the 
blade surface pressure measurements were taken as described in Appendix B of Reference 
10. Inlet total and static pressures were recorded 2 1/3 chord lengths ahead of the blades. 
These pressures were used to non-dimensionalize the blade-surface pressure 


measurements as a coefficient of pressure; 1.e., Cy = (Pjocai - Po) / (Pro - Pwo). 


2; LDV Measurements 


A full listing of the experimental runs is contained in Appendix A. Blade-to-blade 
surveys were conducted at eleven stations ahead of, between, and downstream of the 
blade row. All surveys were taken by collecting 3,000 data points at each location in the 
survey. All LDV data were collected in the “coincidence mode” with a 2.0x10° psec 
window used for validation. Table 1 summarizes the survey stations and locations. Figure 
3 illustrates the same information graphically. The numbers in parentheses on the figure 
indicate the horizontal distance between measurement points. The present study gave an 
average inlet-flow angle of 50.05 degrees. Reference 4 results were for an average B, of 


50.21 degrees. 
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Table 1. LDV Survey Stations and Locations 


Y (inches) 


4d (0.125) 
1b (0.25) 


1 (0.25) 





X (inches) 


Figure 3. LDV Survey Stations Between Blades 6 and 7 


Additionally, eleven specific locations were measured, with an increased number of 
data points, for power-spectrum and auto-correlation analyses. Table 2 below summarizes 
the measurements and the number of data points gathered for each point. Figure 4 
depicts the stations at which these dedicated runs were performed. The points close to the 
blade were chosen to be within the intermittent reverse-flow region, and the outboard 


points were chosen to be outside the reverse-flow region. 


_— ee ee = — ee ———_—_———————————— Se 


Number of Points Acquired 


tation | Point | Location = Location 
(0.0000, -6.2920) 3K, 30K, 50K 
(-1.1331, 4.7905) 3K, 30K 


(0.3654, -3.2920 3K, 30K 
(0.3988, -1.7921) 
(1.0000, -1.7921) 3 


(-0.8394, -4.5420) 3K, 30K 
(0.2500, 0.3620) 


(-0.1225, -3.2920) 3K, 30K 


0K 
(1.5000, 0.3620) 
: 19 23 | (0.2500, 1.0620) 30K 


[19 | 31 | (1.2500, 1.0620) 


Table 2. LDV Power-Spectrum / Auto-Correlation Locations 





The tunnel reference velocity (Vier) for each run was computed as described by 
Elazar in Reference 1. The results are shown in Appendix B. The reference velocity was 
used to non-dimensionalize the data for comparison with previous studies. 

Prior to each run, the probe volume was always referenced to the same location 
between the blades ane an alignment tool. This allowed for a high confidence level in 
specifying station position for comparison with previous studies. The alignment procedure 


was developed and described in detail in Reference 1. 


Y (inches) 





X (inches) 


Figure 4. LDV Power-Spectrum / Auto-Correlation Locations 


a. Inlet Surveys 


The inlet region was surveyed across stations 1, 1b, and 1d (Figure 3). The 
survey spanned three and one-third passage widths in the pitchwise direction (-9 inches to 
+] inch) for station 1, and one and one-third passage widths in the pitchwise direction (-2 
inches to + 2 inches) for stations 1b and 1d. The laser beam system was horizontal when 
measuring at station 1. A 4 degree pitch-up of the LDV table was required for stations 
1b and 1d to avoid blade shadowing of the vertical beam, and to achieve closer access to 
the leading edge. The potential problems due to pitching and yawing the LDV system 
were treated by Hobson and Shreeve in Reference 3. Their analysis showed the maximum 
spatial error due to probe volume orientation to be 0.3 mm, which is the minimum 
diameter of the measurement volume. 

The optics were configured the same throughout the present study with the 
488-nm blue beam measuring the horizontal (U) component and the 514.5-nm green beam 
measuring the vertical (V) component. Frequency shifting was performed as outlined in 


Reference 4 in order to capture the large reversed-flow velocities. 


b. Passage Surveys 

The passage surveys were performed between blades 6 and 7. All passage 
surveys extended 1.8 inches (approximately 60 % passage width) away from the suction- 
side surface of blade 6. The LDV system was pitched up 2 degrees and yawed 4 degrees 
left for stations 2a, 3 and 7 for better access to the suction side of the blade. Stations 10, 


14 and 15 were measured with only 4 degrees of left yaw. 


C. Wake Surveys 


The LDV system was reset to zero pitch and zero yaw for measurements at 
stations 17, 18, and 19. The wake surveys were centered around blades 6 and 7 and 


spanned two passage widths (6 inches). 


d. Laser Power-Spectrum / Auto-Correlation Measurements 


It was desired to expand the analysis of several points throughout the flow 


field by collecting a larger amount of data points. Instead of the standard 3,000 points 


10 


gathered at each location, 30,000 data points were chosen for the locations listed in Table 
2. This choice corresponded to the limiting number of data points that could be acquired 
at positions with low data rates such that the acquisition software reached its maximum 
timeout of 30 minutes. This corresponded to a minimum data rate of 16.7 Hz in order to 
collect 30,000 samples. Stations 2a and 3, near the separation-bubble region, were most 


affected by this limitation. 


3. Flow Visualization 


The laser-sheet flow visualization was photographed through the optical access 
window using a hand-held 8mm camera at 30 frames/second and a VHS camera at 60 
frames/second. The filming was performed at night to reduce glare and optimize the flow 
field image contrast. A total of approximately 15 minutes of footage was obtained of the 
intermittent reversed flow inside the blade passage and of the flow near to the passage 


entrance, including the leading-edge separation bubble. 


1] 





II. COMPUTATIONAL FLUID DYNAMICS 


A. GRID GENERATION 


The grid generation program employed was a modified version of GRAPE (GRids 
about Aerofoils using Poisson’s Equation), [Ref. 11]. This version allowed for the 
generation of a periodic C-type grid for turbomachinery cascades. The namelist input file, 
along with a description of the effort put forth in refining the leading edge geometry for 
the grid, is described in Appendix E. 

The grid, a portion of which is depicted in Figure 5, contained 240 grid points 
wrapped around the blade 1n the streamwise direction and 49 grid points in the blade-to- 
blade direction. The grid was much finer in the normal direction near the surface for 
improved viscous flow resolution, particularly of the leading edge separation bubble. The 
outer, or periodic, boundaries corresponded in size to the cascade wind tunnel as 
described in Experimental Setup. The inlet region was located 30% of an axial chord 
upstream of the leading edge. The wake region extended two chord lengths downstream 


of the trailing edge. 
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Figure 5. Two-Dimensional Periodic C-Grid 


B. COMPUTATIONAL FLOW SOLVER 


Numerical computations were made using Rotor Viscous Code Quasi-3D 
(RVCQ3D) of Reference 8. RVCQ3D solved the thin-layer Navier-Stokes equations by 
employing an explicit, multi-stage, Runge-Kutta algorithm. It could account for the 
effects of rotation, radius change and stream surface thickness variation in blade-to-blade 
turbomachinery flows. The Baldwin-Lomax algebraic eddy-viscosity model [Ref. 12] was 
used to characterize the turbulence within the flow field. 

Computations were performed at the following inlet flow angles: 40, 43, 46, and 
48 degrees. This was done to allow comparisons with existing experimental results. 
Validation of the computed results was achieved by comparing pressure coefficient 
distributions around the blade surface, leading-edge separation bubble sizes, and loss 
computations. 

A typical namelist input file is shown in Appendix F. All solutions were min to 
6,000 iterations and axial-velocity-density ratio (AVDR) was held at 6% for all 
computations. The AVDR was modeled by defining two stream surface points at the inlet 
and exit using the namelist 6 parameter ‘“NMN’ as shown in Appendix F. 

Although the inlet flow angle was specified initially, the flow solver computed to a 
new inlet condition based on the prescnbed exit-static to inlet-stagnation pressure ratio. 
The inlet flow angle varied due to changes, primarily, in the tangential velocity 
component, as the axial velocity did not change due to continuity constraints. 

Solutions were computed using the National Aerodynamic Simulation (NAS) 
Program CRAY C90 supercomputer and the NPS CRAY Y-MP EL98 supercomputer. 
CPU mun time for the NAS C90 was between 3-4 minutes and the NPS EL98 required 
about 55 minutes. Analysis of the solution files was carried out using FAST graphics 
software [Ref. 13] and separate FORTRAN programs given in Appendices H (blade- 


surface pressure distribution) and I (computed inlet-region flow angle). 





A. 


study and Reference 4 is shown in Figure 6. The results, after reconfiguring the tunnel, 
showed reasonable agreement with the pressure distributions obtained earlier at the 
leading edge; particularly on the pressure side. 
trailing edge which was possibly due to a difference in setting of the tunnel’s exit-plane 


tailboards. 


IV. RESULTS AND DISCUSSION 


BLADE-SURFACE PRESSURE COEFFICIENT 


A comparison of pressure coefficient around the blade surface between the present 


distribution. 
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Figure 6. Pressure Distribution Comparison 
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There was some difference near the 


The setting of the tailboards significantly affected the exit static-pressure 


B. LDV SURVEY COMPARISONS 


The data analysis and comparisons which follow were performed without any 
refinement bounds for the velocity histograms. All the results from Reference 4 were re- 


analyzed in this manner. 


1. Inlet Survey (Station 1) 


Figures 7 and 8 show the non-dimensionalized mean velocity and turbulent 
intensities, respectively. The two studies showed variations in the mean velocity of less 
than 2% across three passage widths. The period in the variation in the mean velocity 
corresponded to the 3 inch blade separation, indicating that the blades were felt 30% of an 
axial chord ahead of the leading edge. The U and V turbulence-intensity levels were 
similar, and this indicated the turbulence was isotropic at the inlet. The increase in 
turbulence intensity from 2% to 4% over the extent of the survey is noted. This behavior 


was not observed at lower inlet-flow angles. 


2 Passage Survey (Station 7) 


Mean-flow results for station 7 are presented in Figure 9. There was good 
agreement across the passage which further established repeatability with Reference 4’s 
results. The turbulence intensity in Figure 10 also agreed well and showed anisotropic 
behavior of the turbulence near the blade surface out to 50 % pitch. Beyond this point the 
turbulence became isotropic in the inviscid region, as the levels of turbulence in the U and 
V components were comparable. The data were presented as turbulence intensities, yet 
because the flow field was intermittently reversible the standard deviation of the 
histograms, given by the LDV, did not represent the turbulence in the flow. Rather, these 
should be looked upon as a measure of the unsteadiness in the flow, which was a 


combination of the random turbulence and intermittent reverse flow. 
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Figure 7. Velocity Comparison at Station 1 
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Figure 8. Turbulence Intensity Comparison at Station 1 


19 


Station 7 Mean Velocity 


Non-dimensional Velocity 





X (inches) 


—e— Ganaim U/Vref —8— Ganaim V/Vref —e— Ganaim Utot/Vref 
—e— Present Study U/Vref —a-Present Study V/Vref —e—Present Study Utot/Vref 


Figure 9. Velocity Comparison at Station 7 
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Figure 10. Turbulence Intensity Comparison at Station 7 
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Figures 11 and 12 separate the mean horizontal (U) and mean vertical (V) velocity 
components into their respective positive and negative mean values. The velocity 
histograms for each point across the passage were specially edited such that only negative 
velocities were retained and averaged to create a negative velocity trace and, likewise, 
only positive velocities were retained and averaged at each point to create a positive 
velocity trace. It was apparent from Figures 11 and 12 that the weighted average of the 
negative trace and positive trace would make up the mean velocity for each component 
across the passage as previously shown in Figure 9. 

This analysis showed more of the character of the intermittent reversed flow than 
did a plot of the mean velocity. What occurred in the unsteady stalling process was not 
represented by the arithmetic mean. The envelope about the mean contained the positive 
and negative velocities that truly characterized the flow’s behavior. Since the data were 
taken in coincidence mode, one disadvantage of this form of presentation was that during 
editing of the positive U velocity-component data, all of its negative points were 
discarded, as well as the corresponding V velocity-component data. The opposite but 
equivalent process occurred dunng editing of the V velocity component; such that all the 
data that were left were the velocity components in the first quadrant of a flow that was 
oscillating with components in all four quadrants. It was apparent from the figures that by 
30 % pitch, the flow no longer exhibited reverse velocities in either the U or V 


components. 
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Figure 11. U Velocity Comparison at Station 7 
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Figure 12. V Velocity Comparison at Station 7 
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3: Wake Survey (Station 19) 

The wake surveys were compared at station 19. Figures 13 and 14 show the 
results and indicate good general agreement for the mean velocity and turbulence 
intensity/unsteadiness across two passage widths for the two studies. The periodic 
influence of the blades was strongly evident in the V component of Figure 13. Also, the 
asymmetry of the wake was especially apparent with the larger velocity deficit on the 
suction side of the blade. Figure 14 once again points out the anisotropic behavior of 
turbulence/unsteadiness in this same region. 

Figures 15 and 16 are the results of the special editing as outlined for station 7 in 
Figures 11 and 12. It was of interest to note that the V-component reversed-flow 
velocities at station 19 were of the same magnitude as the V-component reversed-flow 


velocities at station 7. 
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Figure 13. Velocity Comparison at Station 19 
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Figure 14. Turbulence Intensity Comparison at Station 19 
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Figure 15. U Velocity Comparison at Station 19 
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Figure 16. V Velocity Comparison at Station 19 


2 


CG. PASSAGE SURVEY (STATION 3) 


Measurements obtained at 32 locations at station 3 are shown in Figures 17-20. A 
detailed analysis to resolve the unsteady characteristics of the flow was made at the three 
locations (1, 7 and 20) shown in Figure 17 (U component of velocity) and Figure 19 (V 
component of velocity). The histograms of the U and V components of velocity at the 
three locations are shown in Figure 18 and Figure 20, respectively. Point 1 was located 
just off the blade surface. This was inside the leading-edge separation bubble, evidenced 
by the negative mean velocity of all the data points in the histogram (Figures 18a and 20a). 
The second location (point 7) was at 2.5 % pitch. This point was within the intermittent 
reversed-flow region, as indicated by the bi-modal velocity histograms (Figures 18b and 
20b). The third location (point 20) was at 17 % pitch, and was outside the region of 
intermittent reversed-flow. 

Figure 17 shows the mean pitchwise velocity along with the separate averages of 
the positive and negative components. The progression from a negative mean velocity to 
a positive mean velocity was apparent with the crossover distinctly evident in the bi-modal 
histogram at point 7. It was of interest to note that the arithmetic mean for the flow rarely 
occurred. The negative and positive velocity peaks were predominant, which was 
indicative of the unsteady stalled flow. 

The streamwise and pitchwise velocity component behavior at station 3 appear to 
be similar because the blade surface was oriented at approximately 45 degrees to the two 
pairs of laser beams (X-Y axes for the velocity measurements). Figure 21 provides an 
alternate view of the previous data by plotting the pitchwise (U) and axial (V) velocity 


components together. 
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Figure 17. U Velocity Survey at Station 3 
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Figure 18. U Velocity Histograms at Station 3 
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Station 3 V Velocity 
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Figure 19. V Velocity Survey at Station 3 
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Figure 20. V Velocity Histograms at Station 3 
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Figure 21. U- V Velocity Correlation Plot 
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D. POWER-SPECTRUM / AUTO-CORRELATION ANALYSIS 


Power spectral density and auto-correlation analyses were performed on the data 
collected as previously indicated in Table 2. Only station 3 (point 7) and station 7 (points 
3 and 20) have been included in Appendix C. Therein are contained the velocity 
histograms and time-history plots for these points. This analysis was performed using 
TSI’s FIND software, version 4.03. The settings chosen for the analysis were the 
following: Window Type: Hanning, Fast Fourier Transform (FFT) Method: Random 
Correlation, and Correlation Type: Random Correlation. In all cases the default lagged- 
time slot (At) was used 1n order to satisfy the Nyquist sampling criterion. 

It was initially believed that a large data set was required to properly carry out this 
type of analysis. The data sets acquired for station 3, point 7, included 3K and 30K data 
points. The comparison of the two sets showed similarity in their behavior. The same can 
be said for station 7, points 3 and 20, which were compared between 3K, 30K, and 
Ganaim’s 3K data points. The only point of recurring significance throughout the power 
spectrum and auto-correlation plots was in the region between 20-25 Hz as shown in 
Figure 22 (30K points, present study), and Figure 23 (3K points, Reference 4) for station 
7, point 3. No further interpretation of the results obtained was attempted due to the many 
limitations of discretely-sampled, randomly-arnving data points. Additionally, since this 
routine involved averaging the data over the lagged-time slot (At), there was some 


question whether the true character of the flow had been captured. 
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Figure 22. Power Spectrum for 30K Points 
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Figure 23. Power Spectrum for 3K Points 
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The use of larger data sets did not appear to be of significance in the final analysis. 
The primary limitation turned out to be the data sampling rate. The maximum frequency 
for spectral analysis was limited to 1/2 the data rate due to the Nyquist criterion, which 


prevents aliasing of the data. Table 3 shows the data rates and maximum frequencies for 


the spectral analysis plots of Appendix C. 
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Table 3. Spectral-Analysis Maximum Frequencies 


E. FLOW VISUALIZATION 

Visualization of reversed-flow velocities within a blade passage was documented 
on VHS videotape. The tunnel was run at a plenum chamber pressure of 12 inches H20 to 
match the conditions during the LDV measurements. The shutter speeds for the 8mm 
recording.were varied from 1/250th to 1/1000th and 1/4000th of a second with the rate of 
30 frames/second remaining constant. It was found that the slower shutter speed of 
1/250th of a second was better able to capture the streaklines of the intermittent reversed- 
flow due to longer particle traces for each image. 

Upon replaying the videos at normal and slower speeds, several observations were 
made for both the present study and the flow visualization performed in Reference 4. The 
reversed-flow region appeared to be confined to the blade passage and did not extend 


upstream of the leading edge. The leading-edge separation bubble trapped and 
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accumulated seed particles, due to injection from the reversed flow. This seed within the 
bubble was able to sustain itself even after the fog was turned off. The bubble also 
appeared to pulsate in size in accordance with the intermittent reversed-flow. There was 
some three-dimensionality to the reversed flow evidenced by small spanwise velocity 
traces. Since the laser sheet was only projected between two blades, it was not possible to 
determine whether the cascade was experiencing intermittent stall from blade passage to 


blade passage, or if the entire cascade was stalling simultaneously. 
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F. COMPUTATIONAL RESULTS 


1. Comparison of Blade-Surface Pressure Distributions 


The predicted and experimental blade-surface pressure distributions are shown for 
40, 43, 46, and 48 degrees in Figures 24-27 respectively. Experimental data for the 40 
and 43 degree case were from Reference 14. The 46 degree case used Reference 15 
experimental data and the 48 degree case used Reference 3 data. The AVDR parameter 
was held constant at 6% for all solutions, which showed good agreement for the first three 
cases. The 48 degree prediction over-estimated the experimental values aft of 60% chord. 
This was most probably due to under-estimating the stream-tube thickness variation at this 
inlet flow angle. Increasing the AVDR, as inlet flow angle increased, improved the 
agreement of the computational results; however, this was not the case experimentally, as 
documented in References 14, 16, and 17. 

The computed static-pressure recovery far downstream of the trailing edge and 
that measured (from Reference 18) showed good agreement for the 40 and 43 degree 
cases. The computed value for 40 degrees was 0.313 as compared with a measured range 
of 0.32 - 0.34. The computed value for 43 degrees was 0.357 as compared with a 
measured range of 0.37 - 0.39. The computed static-pressure rise far downstream is 
relevant to accuracy in predicting the losses using viscous codes, and it was shown here to 
agree well with experiments near design incidence. 

For B; = 48 degrees, the predicted flow field pressure contours between the blade 
rows are shown in Figure 28. The pressure field was continuous at the grid boundaries 
showing good periodicity; however, there was some indication of grid contamination of 
the pressure field at the trailing edge cut. The suction peak and stagnation point were 
both clearly evident. The static-pressure nse, or adverse-pressure gradient, towards the 
trailing edge on the suction side was also apparent. Investigation of the velocity field in 
this region showed that there was flow separation aft of mid-chord, followed by re- 


attachment prior to the blunt trailing edge. Similar behavior occurred for the 46 degree 
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solution. Experimentally this was not the case, as the flow remained attached on the aft 
section of the suction surface. This indicated that the code may be predicting incipient 
stall at 48 degrees (or at a lower inlet-flow angle of 46 degrees); whereas, experimentally 
this was only found at 50 degrees. It appeared that an unsteady solution may have been 
computed as evidenced by the climbing density residuals, which are shown in Appendix G, 
and the coalescence of the pressure contour lines throughout the flow field (Figure 28). 


The 40 and 43 degree solutions did not show any suction-side separation on the aft 


portion of the blade, in agreement with experiments. 
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Figure 24, 40 Degree Pressure Distribution 
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Figure 25. 43 Degree Pressure Distribution 
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Figure 26. 46 Degree Pressure Distribution 
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Figure 27. 48 Degree Pressure Distribution 
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Ds Leading-Edge Separation Bubble 


A leading-edge separation bubble was predicted for all computational solutions. 
The location of the computed separation point was insensitive to varying inlet flow angle. 
The re-attachment point for each case was determined graphically. Figure 29 shows the 
bubble region for B; = 48 degrees in which the recirculating flow was evident. A 
comparison of computed and experimental re-attachment points is shown in Figure 30. 
The experimental data was from Reference 19, obtained using the china-clay technique. 
The plot showed an under-prediction of greater than 50% for the re-attachment locations, 
which is typical of bubble computations. The Baldwin-Lomax turbulence model does not 
account for freestream turbulence, which, as pointed out in Reference 3, is necessary in 
order to properly characterize the development of the blade surface boundary layer, or the 
occurrence of separation. The two most likely reasons for the shortcoming of the 
computations are firstly, the limitations of the Baldwin-Lomax turbulence model and 
secondly, the lack of grid refinement in the streamwise direction aft of the leading edge. 
The trend of increasing downstream distance for re-attachment, as inlet flow angle 


increased, appeared to be properly computed. 
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Figure 29. Leading-Edge Velocity Plot 
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Figure 30. Separation-Bubble Reattachment 
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3. Losses 

The final assessment of RVCQ3D was its loss-prediction capability. Figure 31 
plots measured losses as reported in References 14, 16 and 17 with those computed in the 
present study. The first two computed cases of 40 and 43 degrees were within 5% of the 
measured results. The 46 degree case showed a sharp departure from the measured 
results, falling short by 45%. An attempt to increase the AVDR parameter for this case 
led to a small increase in the loss prediction but with a marked degradation of the blade 
pressure distribution. 

The output file for each solution contained a listing of the Baldwin-Lomax 
turbulence-model behavior by printing values of turbulent eddy viscosity for each grid 
point. Upon interrogation of this output file, the behavior of the Baldwin-Lomax 
turbulence model was seen to exhibit certain characteristics. Firstly, for all cases the 
turbulence model turned on inside the separation bubble region. Secondly, for the 40 
degree case, the model did not turn on until two grid points inside the bubble region, 
indicating a laminar region in the forward section of the bubble. The model then turned 
off again at the re-attachment point and turned back on further downstream. Thirdly, for 
48 degrees the model turned on at the separation point and remained on past the re- 
attachment point to the trailing edge. The transition criteria used for all the computations 
was that once the turbulence model computed a turbulent eddy viscosity (Lt7) which was 
14 times larger than the laminar viscosity (iL), this value of pr was used to compute an 
effective viscosity (le = Ui + Hr ). Otherwise the effective viscosity was simply the 


laminar viscosity. 
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Figure 31. Loss Comparison 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Experimental flow field measurements near stall were successfully repeated for a 
controlled-diffusion compressor cascade using a two-dimensional laser-Doppler 
velocimeter system. The nominal inlet flow angle achieved was 50.05 degrees at a 
constant Mach number of 0.25 and Reynolds number of 700,000. The measurements 
were useful in that they showed repeatability with previous work by Ganaim, and 
validated a test case for two-dimensional unsteady separated flows. 

Analysis of the data was presented through a unique editing process which 
captured the character of the intermittent reversed flow. This process more properly 
represented the stalling region as a combination of primarily positive and negative flow, 
since the time-mean velocity rarely occurred in the intermittent region. 

Power-spectrum and auto-correlation calculations were completed on the data for 
eleven locations. Interpretation of these results was limited to their general behavior due 
to averaging of the data by the software package. Three locations were chosen for 
analysis and the results were found to depend greatly on the data acquisition rate. When 
expecting to perform power-spectrum analysis on experimental data, it 1s important to 
ensure that the data rate is high by repeated cleaning and proper re-fitting of the optical- 
access window. 

Numerical computations showed good predictability of the blade pressure 
distribution up to 46 degrees (approximately 6 degrees above design incidence). The 48- 
degree case began to deviate by up to 50% from measured values aft of two-thirds chord. 
A case-by-case basis for adjusting AVDR may be necessary to bring the computed results 
into better agreement with the measurements at higher inlet flow angles. 

A leading-edge separation bubble was formed in all computations in accordance 
with experiments. The extent of the computed bubble fell short of that seen in 
measurements. Computed losses agreed well at low incidence angles. As the inlet flow 
angle increased above 46 degrees, predicted losses deviated by 45% from the 
measurements. The use of a constant AVDR for numerical solutions greatly simplified the 


process. This was found to be acceptable for analyzing near on-design inlet flow angles. 
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Future computations of the flow through the cascade with a quasi-3d code should 
focus on achieving the following: 

1. Continue numerical solutions to determine the optimum AVDR as inlet-flow 
angle varies, and compare with those measured in the cascade. 

2. Quantify what effects adjusting AVDR to match pressure distributions has on 
the computed losses. 

3. Increase the grid size in both the I and J directions for better characterization of 
the flow near the blade surface and around the leading edge separation bubble. 

4. Replace the Baldwin-Lomax turbulence model with a two-equation turbulence 
model to account for freestream turbulence intensity. 

5. Perform time-accurate solutions at higher incidence angles to see if 


unsteadiness is being properly modeled. 
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APPENDIX A: EXPERIMENTAL-MEASUREMENTS LOG 


A detailed log of the atmospheric conditions, tunnel plenum chamber 


conditions and LDV traverse-mechanism geometry for each aan asa run. 
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APPENDIX B: CALCULATION OF INLET-FLOW REFERENCE VELOCITY 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


OUTPUT FROM PROGRAM CALIBRATE 
cccccccccceccccccecccccccccececccccccccccccccccccccccccccccceccccccc 


LEAST SQUARES STRAIGHT LINE CURVE FIT IS USED 
TO DETERMINE TUNNEL CHARACTERISTICS AT DIFFERENT SPEEDS 


NEWTON S METHOD IS USED TO DETERMINE THE REFERENCE VELOCITY 
FROM THE RECORDED AMBIENT PRESSURE AND TUNNEL PLENUM 
PRESSURE AND TEMPERATURE 


BEGIN DETERMINING TUNNEL CHARACTERISTICS 
FROM THE FOLLOWING MEASURED VALUES 


AXIAL VEL. TANGENTIAL VEL. AMBIENT PRESS. PLENUM PRESS. PLENUM TEMP. 
M PER SEC. MPER SEC. INCHES MERCURY INCHES WATER _ DEG. C. 

18.3100 22.1350 29.9702 1.7000 13.8888 

32.4070 38.5280 29.9702 4.7000 13.8888 

42.5440 90.1820 29.9702 8.1000 15.5595 

47.2920 96.3530 29.9702 10.0000 16.6666 

51.7870 61.9370 29.9702 12.1000 ezeee 

959.3660 66.7010 29.9702 14.0000 18.3333 


CALCULATED VALUES FOR THE TUNNEL CONFIGURATION 


TOTAL VELOCITY 


0.28726543910E+02 
0.50345013984E+02 


MACH NUMBER 


0.37818777200E-01 
0.66279705382E-01 


MACH NUMBER FUNCT. PRESSURE RATIO 


0.49880294974E-02 
0.15207192279E-01 


0.43443528269E-02 
0.11919477085E-01 


0.65789247298E+02 0.86361830486E-01 
0.73567614295E+02 0.96387235172E-01 
0.80734660079E+02 
0.86685 739064E+02 


0.25620261319E-01 
0.31766759119E-01 
0.10567615343E+00 0.38003966212E-01 
0.11324925118E+00 0.43463396462E-01 


0.20366465634E-01 
0.25024235973E-01 
0.30121036781E-01 
0.34686726827E-01 


CALLING LEAST SQUARES SUBROUTINE 
TO DETERMINE THE PRESSURE RATIO AS A FUNCTION OF MACH NO. PARAM 


PRESSURE RATIO = A1 * ANUX + AO 
MATRIX EQUATION 


0.60E+01 0.16E+00A0 0.13E+00 
0.16E+00 0.53E-02A1 0.42E-02 


A1= 0.78899597626E+00 AOQ= 0.16213280706E-03 
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REFERENCE CONDITIONS FOR EACH RUN 


1 


PRESSURE RATIO = 


AMB PRESS PLENUMPRESS PLENUM TEMP RUN NAME 

INCHES Hg INCHES H20 DEGREES C 
29.9295 11.9000 18.3333 1114S51_A.PRN 
30.1331 12.0000 14.4444 ST1_3K.PRN 
30.1331 12.0500 14.4444 ST1_30K.PRN 
30.1331 12.0000 15.5555 ST1_SOK.PRN 
30.1127 12.0500 16.6666 1121S_1B.PRN 
30.1127 12.1000 18.3333 1121S_1D.PRN 
30.2349 11.9500 16.9444 1122S_2A.PRN 
30.1331 12.0000 VET CT? ST2A_3K.PRN 
29.9295 12.1500 14.4444 ST2A_30K.PRN 
30.1942 11.8500 17.2222 1123S_3.PRN 
30.1331 12.0000 16.6666 ST37_3K.PRN 
29.9295 12.0000 16.6666 ST37_30K.PRN 
30.2145 12.0000 20.0000 1123S_7.PRN 
30.1331 12.1000 15.5555 ST73_3K.PRN 
29.9295 12.0000 17.2222 ST73_30K.PRN 
30.1331 12.1000 14.4444 ST7_3K.PRN 
30.0313 11.9000 14.4444 ST7_30K.PRN 
30.0517 12.1000 17.2222 ST10_EDG.PRN 
30.0313 12.0500 16.6666 ST10_OUT.PRN 
30.1738 12.1000 15.5555 1127S _14.PRN 
30.2145 12.0500 16.1111 1128S_15.PRN 
30.2145 12.1000 17.2222 1128S_17.PRN 
30.0313 11.9000 16.1111 ST17_IN.PRN 
30.0313 11.9000 16.6666 ST17_OUT.PRN 
30.2349 11.9000 18.3333 1128S_18.PRN 
30.2145 12.0000 19.4444 1128S_19.PRN 
30.0517 12.2500 16.6666 ST19_IN. PRN 
30.0517 12.1000 16.6666 ST19_OUT.PRN 


0.02968 MACH NUMBER PARAMETER = 0.3741E-01 


RUN NAME = 111481_A.PRN 
BEGIN NEWTON ITERATION 


ITERATION NUMBER 1 MACH NO. PARAM. = 0.105475 ERROR TERM = 0.6522E-03 
ITERATION NUMBER 2 MACH NO. PARAM. = 0.104822 ERROR TERM= 0.1641E-05 
ITERATION NUMBER 3 MACH NO. PARAM. = 0.104821 ERROR TERM = -0.3698E-09 
ITERATION NUMBER 4 MACH NO. PARAM. = 0.104821 ERROR TERM= 0.8591E-13 


VREF = 80.23415022742 


l= 28 
PRESSURE RATIO = 0.03004 MACH NUMBER PARAMETER = 0.3787E-01 
RUN NAME = ST19_OUT.PRN 


BEGIN NEWTON ITERATION 
ITERATION NUMBER 1 MACH NO. PARAM. = 0.105777 ERROR TERM= 0.2921E-03 
ITERATION NUMBER 2 MACH NO. PARAM. = 0.105485 ERROR TERM= 0.2872E-06 
ITERATION NUMBER 3 MACH NO. PARAM. = 0.105485 ERROR TERM = -0.6810E-10 
VREF = 80.51152210108 


54 


EXPERIMENT REFERENCE 
VELOCITY 


NUMBER 


OcoOnNO  & WADA = 


80.2342 
79.7607 
79.9273 
79.9146 
80.2627 
80.6606 
79.8038 
80.2216 
80.5328 
79.5598 
80.0682 
80.3411 
80.4184 
80.2481 
80.4180 
80.0935 
79.5615 
80.5887 
80.3718 
80.1937 
80.0501 
80.3705 
79.7917 
79.8682 
79.8265 
80.3422 
81.0107 
80.5115 


FILE NAME 


1114s1_A.PRN 
ST1_3K.PRN 
ST1_30K.PRN 
ST1_50K.PRN 
1121S_1B.PRN 
1121S_1D.PRN 
1122S _2A.PRN 
ST2A_3K.PRN 
ST2A_30K.PRN 
1123S_3.PRN 
ST37_3K.PRN 
ST37_30K.PRN 
1123S_7.PRN 
ST73_3K.PRN 
ST73_30K.PRN 
ST7_3K.PRN 
ST7_30K.PRN 
ST10_EDG.PRN 
ST10_ OUT.PRN 
1127S_14.PRN 
1128S_15.PRN 
1128S_17.PRN 
ST17_IN.PRN 
ST17_OUT.PRN 
1128S_18.PRN 
1128S_19.PRN 
ST19_IN.PRN 
ST19_OUT.PRN 
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APPENDIX C: VELOCITY HISTOGRAMS, TIME-HISTORY PLOTS, POWER- 


SPECTRUM PLOTS, AND AUTO-CORRELATION PLOTS 





1. Station 3, Point 7 
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b. 30K Data Points 
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2: Station 7, Point 3 
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APPENDIX D: REDUCED DATA FROM STATIONS 1, 3, 7, 19 


X Gn) | Y Gn) | U/ Vref [ Ve Vref [UV Mag/ Vref | U- Turbulence | V- Turbulence | Reynolds | Correlation 
meee ? Stress Coefficient 


| 1.00 | 629 | 0.7560 0.6273 0.9824 4.5167 4.2288 0.5645 0.0459 
0.7527 0.6348 0.9847 4.2846 4.1545 0.1621 0.014] 
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Pitchwise Survey at Station 19 
[CdSe Flow Average] —=«d;_S*d;sCAn erage of Negative | Average of Postve 


Vref 


59995 0.6604 077056 
38750 0.6619 
57499 | 1.0620 | 00195 | 0sae1 | 05864 | 85960 | 256453 | 00069 | 00846| o.0se4 | 0.6254 
3.6250 | 1.0620 | 00197 | 0.5217 | 05220 | 9.1367 | 245519 | 00600 |-o0as1| 0.0751_| 05499 
55001 0.4960 
53751 0077 
| 03312 
3.1251 
5 
2.8750] 10620 [0.055 0.7741 | 0.764 | 63753 | 110721 | 0.0000 | 0.0000 | a0748 | 0.7823 
27498 | 1.0620 | 0.0458 | 0.8239| 0.4251 | 5.0207 | 5.1598 | 0.0000 | 0.0000 | 0.0367 | 08354 | 
08273 
0361 
078268 


© 


1.8750_| 1.0620 0.8332 
[7500 | 1.0620 006 08173 
6250__| 1.0619] 00368 | 08066 | 08075 | 3.1928 | 10.1344 | 0.0000 | 0.0000 0.8132 
13000 | 1.0620 Doasi| 0.0707 | 0.7438 
13750 | 1.0620 0.0590] 0.0711 | 0.7229 
0-736 

; 07198 
; : : 0 


2 08273 


oO; o 
© 


S 


E 


© 


.00 0.6013 
0.5917 


: 0 
0.7499 | 1.0620 0.6326 
0.6250 | 1.0620 73.6395 0.5647 
0.5000 | 1.0620 04977 
03750 | 1.0620 [00321 0.4085 
02500 | 1.0620 | 0.0439 | 03235 | 03265 | 12.6135 | 169481 | 00607 |-0.0682| 0.1039 | 03344 
0.1250 0.3390 
2.0001 0.5487 
2.1230 070600 0.8025 
22501 0.8430 
D375 
2.4999 
2.6250 
0.7500 
0.8750 
“7.0000 
1251 
712500 
13750 
“5000 
6250 08119 
“1.7500 0.0646 | 0.7614 
“18750 | 1.0620 [0037907059 | 07069 | Tass | 223243 | 0s |oorm| 0007 | 078i 
2.0000 | 1.0620 [oosre [asso] 0.46927 | 76868_| 237145 | 0.0669 |-o.o7e8] 0.0598 | 0.7296 


©& —_— — 
00 r — fa a 
31815 8 8 
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APPENDIX E: CFD GRID GENERATION 


The FORTRAN grid generation program employed was the modified version of 
GRAPE (GRids About Aerofoils using Poisson’s Equation) by R. V. Chima, [Ref. 11]. 
Prior to this study, the grid’s leading edge was blunt or flat-faced, which led to 


inaccuracies in the computed solutions. It was determined that GRAPE properly 


smoothed the points around the leading edge, and that it required definition of the leading 


edge with more user-supplied coordinates. Careful attention was given to properly 


blending the top and bottom grid surfaces into the 0.045 inch radius circular arc that 
defines the leading edge. This was done by ensuring that first derivative continuity was 
within a few percent difference. The final coordinates for the CD blade definition are 


contained below in namelist &GRID3. 


&GRID1 

JMAX=240,KMAX=49, NTETYP=3,NAIRF=5, JAIRF=163,NIBDST=7, 

DSI=0.0002, JTEBOT=60,JTETOP=181,XLE=0.0,NOBSHP=7,XTE=1.0, _ 

ALAMF=0.0,ALAMR=0.0,XLEFT=-.3,XRIGHT=3., 7 

RCORN=.10,NOUT=4,NORDA=5,3,MAXITA=100, 100, 

&END 

&GRID2 

AAAI=.3, BBBI=.3, DSOBI=0.01,ROTANG=-14.4, 

XTFRAC=1 .5,PITCH=.5988, DSRA=.4920, DSLE=.0008, DSTE=.001, 
=28, NTE=18, WAKEP=1.,OMEGR=1.0,OMEGS=1 .0,OMEGP=1.0,OMEGQ=1.0, 

&END 

&GRID3 AIRFX= 

0.9996, 0.9971, 0.9929, 0.9876, 0.9830, 0.9803, 0.9776, 0.9749, 


0.9601, 0.9453, 
0.8419, 0.8271, 
0.7238, 0.7090, 
0.6056, 0.5908, 
0.4874, 0.4726, 
0.3693, 0.3545, 
0.2511, 0.2364, 
0.1329, 0.1181, 
0.0244, 0.0114, 
0.0033, 0.0021, 
0.0010, 0.0019, 
0.0591, 0.0738, 
0.1772, 0.1920, 
0.2954, 0.3102, 
0.4136, 0.4283, 
0.5317, 0.5465, 


0.9305, 0.9157, 0.9010, 0.8862, 
0.8124, 0.7976, 0.7828, 0.7681, 


0.6942, 0.6795, 
0.5760, 0.5612, 
0.4579, 0.4431, 


0.3397, 0.3250, 0.3102, 0.2954, 


0.2216, 0.2068, 
0.1034, 0.0886, 
0.0097, 0.0080, 
0.0011, 0.0005, 
0.0030, 0.0044, 
0.0886, 0.1034, 
0.2068, 0.2216, 
0.3250, 0.3397, 
0.4431, 0.4579, 
0.5612, 0.5760, 


0.6647, 
0.5465, 
0.4283, 


0.1920, 
0.0738, 
0.0063, 
0.0001, 
0.0045, 
0.1181, 
0.2364, 
0.3545, 
0.4726, 
0.5908, 


0.8714, 
0.7533, 
0.6352, 
0.5169, 
0.3988, 
0.2807, 
0.1772, 0.1624, 
0.0591, 0.0443, 
0.0047, 
0.0001, 
0.0177, 
0.1329, 
0.2511, 
0.3693, 
0.4874, 
0.6056, 


0.6499, 
0.5317, 
0.4136, 


0.0004, 
0.0310, 
0.1477, 
0.2659, 
0.3840, 
0.5022, 
0.6204, 


ws 


0.8567, 
0.7385, 
0.6204, 
0.5022, 
0.3840, 
0:2659, 
0.1477, 
0.0343, 


0.0443, 
0.1624, 
0.2807, 
0.3988, 
0.5169, 
0.6352, 


0.6499, 
0.7681, 


0.6647, 
0.7828, 
0.9010, 


0.6795, 0.6942, 
0.7976, 0.8124, 
0.9157, 0.9305, 
0.9956, 0.9988, 


0.0012, 
0.0178, 0.0218, 
0.0431, 0.0453, 
0.0575, 0.0589, 


0.0661, 
0.0714, 
0.0688, 
0.0509, 
0.0238, 
0.0001, 
0.0047, 
0.0157, 
0.0613, 
0.1074, 
0.1345, 
0.1413, 
0.1320, 
0.1090, 
0.0784, 
0.0456, 
0.0219, 


0.7090, 
0.8271, 
0.9453, 
1.0000, 


0.0256, 
0.0474, 
0.0601, 


0.7238, 
0.8419, 
0.9601, 
0.9996, 


0.0000, 0.0008, 0.0018, 


0.0291, 
0.0494, 
0.0612, 


0.0669, 0.0676, 0.0683, 


0.0717, 
0.0675, 
0.0478, 


0.0719, 
0.0659, 
0.0446, 


0.0719, 
0.0641, 
0.0413, 


0.0202, 0.0165, 0.0127, 


0.0001, 
0.0063, 
0.0167, 
0.0680, 
0.1120, 
0.1363, 
0.1411, 
0.1297, 
0.1054, 
0.0744, 


0.0005, 
0.0080, 
0.0168, 
0.0743, 
0.1163, 
0.1379, 
0.1406, 
0.1272, 
0.1018, 
0.0704, 


0.0011, 
0.0097, 
0.0244, 
0.0804, 
0.1204, 
0.1391, 
0.1399, 
0.1246, 
0.0981, 
0.0663, 


0.0415, 9.0373, 0.0331, 
0.0176, 0.0124, 0.0093, 


28 


0.7385, 0.7533, 
0.8567, 0.8714, 
0.9749, 0.9802, 


0.0028, 0.0038, 
0.0323, 0.0353, 
0.0513, 0.0530, 
0.0623, 0.0633, 
0.0690, 0.0697, 
0.0717, 0.0713, 
0.0621, 0.0597, 
0.0379, 0.0344, 
0.0088, 0.0061, 


0.0114, 

0.0319, 0.0391, 
0.0863, 0.0920, 
0.1240, 0.1272, 
0.1401, 0.1407, 
0.1389, 0.1375, 
0.1218, 0.1188, 
0.0943, 0.0904, 
0.0622, 0.0580; 
0.0289, 0.0274, 


APPENDIX F: RVCQ3D CODE INPUT 


A sample namelist input file for RVCQ3D is shown below. Namelist 6 
contains the definition for a linear 6% streamtube contraction which accounts for AVDR 
and the “Quasi-3D” effect in the code. All solutions were run to 6,000 iterations and 


consumed, on the average, only 55 minutes on the NPS CRAY (“Sinus”) and 4 minutes 


on the NAS CRAY (“Von Neumann’). 


| ‘SANGER CONTROLLED DIFFUSION BLADE' 
&NL1 M=240,N=49,MTL=60,MIL=108 &end 


&NL2 NSTG=+4,IVDT=1 IRS=1,EPI=.30,EPJ=.40,CFL=4.0 AV2=0., 
AV4=1. & END 


&NL3 IBCIN=1,IBCEX=1 ITMAX=6000, TIRES TI=0,IRESTO=1,IRES=1 ICRNT=10000, 
IXRM=0 &END 


&NL4 AMLE=0.24,ALLE=43.0, BETE=0,PRAT=0.9750,POIN=1.0, TOIN=1.0;G=1.40, 
&END 


&NLS ILT=2,JEDGE=30, RENR=1.0e6,PRNR=.70,PRTR=0.9,TW=0.0, VISPWR=.666666, 
CMUTM=14.0 &END 


&NL6 OMEGA=0.0,NBLADE=1,NMN=2 &end 
| ~,3000 3.000 

1.0000 1.000 

1.0000 .9434 





94 


APPENDIX G: DENSITY RESIDUALS FOR CFD SOLUTIONS 


40 Degree Solulion (6% AVOR) 
88-05 


*AU22.1e3° —— 


8a-05 }- 
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hii ee ID 
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3000 4000 


0 1000 2000 
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AB Degree Solullon (6% AVOR) 
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{000 2000 3000 dooo 5000 OVO 
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APPENDIX H: Cp VERSUS X/C FORTRAN PLOTTING PROGRAM 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C Program To Read The Output From Rvcq3d.F And Grape.F C 
C And Generate A Data File For A P/P0 Vs Chord Plot C 
C Around The Airfoil. Aug.22, 1991 C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C IMPLICIT REAL*8(A-H,O-Z) 
DIMENSION Q(300,300,4),X(300,300), Y (300,300) 
DIMENSION U(300,300), V(300,300),P(300, 300) 


C CALCULATE THE CRITICAL VELOCITY 
cepe=1.00 
rgas=1.0 
pi=1.0 
ti=1.0 
PRAT=0.685 
rhoi=pi/(rgas*t)) 
ceve=cepe-rgas 
g=1.40 
gp=gtl 
cstar=sqrt(2*g*pi/(gp*rhoi)) 


ISTART=60 
IFINIT=1381 


READ(7,*)NI,NJ 

print *,ni,nj 

READ(7,*)((X(,J),IF1,ND,J=1,NJ), 
& ((Y(LJ,I=1,ND,JH1,NJ) 

READ(7,*)MTL,MIL 


READ(3,*)NI,NJ 
READ(3,*)FSMACH,ALF,RE,TIME 
READ(3,*)(((Q(L,J,K),I=1,ND,J=1,NJ),K=1,4) 
GAMMA=1.4 
DO 1234 I=1,NI 
DO 4321 J=1,NJ 
U_)=Od,J,2VQ(LJ,1) 
VIJ=Q.,J5,3 VOCS, 1) 
P(L,J)=(GAMMA-1)*(Q(L,J,4)-.5*Q(LJ,1)*(U,J)**2+V(1,J)**2)) 
&  *RHOI*CSTAR**2 
4321 CONTINUE 
1234 CONTINUE 

C 

C DETERMINE THE XMIN AND XMAX GRID POINT POSITION 
XMIN=0.0 
XMAX=0.0 
DO 2 I=ISTART,IFINIT 
IF(X(1,1).LE.XMIN)THEN 

=| 
PRINT *,IMIN 
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XMIN=X(L 1) 
END IF 
IF(X(1,1).GE.XMAX)THEN 
IMAX=I 
c PRINT *,IMAX 
XMAX=X(L1) 
END IF 
2 CONTINUE 
C 
C FREE STREAM STATIC PRESSURE USES THE XMIN GRID POINT 
PO=P(IMIN, NJ)+.5*Q(IMIN, NJ, 1)*(U(IMIN, NJ)**2+V(IMIN,NJ)**2) 
& *RHOI*CSTAR**2 


i=imin 
Jn) 
tO=1/rgas*((g-1)*q(ij,4)/q@y, 1)-(g-1)**2/(2*g)* 
& (q(ij,2)**2+q(ij,3)**2)/qaJ,1)**2) 
t0=t0*cstar* *2 
print *,p0,t0,p(imin,nj),q(imin, nj, 1),u(imin,nj),v(imin,nj) 


CALCULATE THE MASS FLOW RATE FROM THE INPUT DATA 
NORMALIZED WITH THE INLET AREA 
dmass=pi/(rgas*ti)*(1/prat)* *(-1/g)*sqrt(g*(g-1)*rgas*ti/2* 
&  (1-1/(1/prat)**((g-1)/g))) 
print *,'cal. inlet mass flow rate =',dmass,'* area’ 


PS Le 25) 


qQ © 


dmasse=q(imin, nj, 1)*RHOI*u(imin, nj) 
print *,’comp. inlet mass flow rate / area =',dmasse 


Oy @re 


PRINT OUT DOWN STREAM CONDITIONS 
pl=p(1,1)+.5*q(1,1,1)*(ud1, 1)**2+v(1, 1)**2)*RHOI*CSTAR**2 
t1=1/rgas*((g-1)*q(1,1,4)/q(1, 1, 1)-(g-1)**2/(2*g)* 

& (q(1,1,2)**2+q(1,1,3)**2)/q(1,1,1)**2) 
tl=tl1*CSTAR**2 

print *,'downstream condition’ 
print *,p1,t1,p(1,1),q(1,1,1),ud,1),vq, 1) 


Ngee a 


CALCULATE THE MASS FLOW RATE AT EXIT 
smass=0.0 
ra=0.0 
ba=0.0 
do 21 j=1,nj-1 
smass=smass+(q(1,j,1)+q(1,j+1,1))*(u(1 y)+ud1,J+1))* 
& RHOI*CSTAR*(abs(y(1,j+1)-y(1,j)))*0.25 
ra=rat+(q(1,j,1)+q(1j+1,1))*(abs(y(1j+1)-y(1,j)))*0.5* 
& RHO! 
ba=batabs(y(1.jJ+1)-y(1,j)) 
21 continue 
do 22 j=1,nj-1 
smass=smass+(q(ni,j,1)+q(nij+1,1))*(u(ni,j)+u(ni,J+1))* 
& RHOI*CSTAR*(abs(y(ni,j+1)-y(n1,j)))*0.25 
ra=rat(q(nl1j,1)+q(mij+1, 1))*(abs(y(ni,j+1)-y(ni,j)))*0.5 
& *RHOI 
ba=batabs(y(ni,j+1)-y(ni,j)) 
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22 continue 
C 
C AVERAGE VELOCITY 
va=smass/ra 
dr=ra/ba 
print *,'vel. downsteam,vd =',va 
print *,'ave. density,dr =",dr 


CALCULATE THE INPUT INLET VELOCITY 
NORMALIZED BY CSTAR 
cO=sqrt(g*rgas*ti) 
vin=sqrt(g*(g-1)*rgas*ti/2*(1-(1/prat)**((1-g)/g)))/c0 
print *,'input inlet velocity / cO =',vin 


C1'O.@: 2) 2 


OCG 


CALCULATE THE CORRECTED DOWNSTREAM VELOCITY 
do 23 j=nj,1,-1 
vac=u(ni,j)*vin/va 
alpha=atan(abs(v(ni,j)/u(ni,j))) 
vic=vac*tan(alpha) 
vc=sqrt(vac* *2+vtc**2) 
yy=abs(y(ni,nj)-y(1,nj)) 
‘s write(24,*)abs(y(ni,j)-y(ni,nj))/yy,ve 
23 continue 
G 
do 25 j=2,nj 
vac=u(1,j)*vin/va 
alpha=atan(abs(v(1,j)/u(1J))) 
vtc=vac*tan(alpha) 
ve=sqrt(vac**2+vic**2) 
write(24,*)abs(y(1,j)-y(ni,nj))/yy,ve 
25 continue 


Q 


CHORD=ABS(X(IMAX, 1)-X(IMIN, 1)) 


ODO. 


DO 1 I=ISTART,IFINIT 
DIST=ABS(X(L, 1)-X(IMIN,1)) 
XS=DIST/CHORD 
CP=1.0+(-P0+P(I,1))/(.5*Q(IMIN,NJ, 1)*(U(IMIN,NJ)**2 
& +V(IMIN,NJ)**2)*RHOI*CSTAR**2) 
vver=sqrt(u(i, 1)**2+v(i, 1)**2) 
C THE SURFACE POINT P(,1) ALWAYS EQUAL 0? 
PT=P(1,1)/PO 
c  WRITE(63,*)XS,VVCR 
c WRITE(64,*)XS,PT 
WRITE(65,*)XS,CP 
1 CONTINUE 
C 
STOP 
END 
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APPENDIX I: UPSTREAM FLOW ANGLE FORTRAN PROGRAM “angle.f” 


ececeeceeceeeceeeceececececeeecceecececeeeeeceeceeeceeececceceeeceeececececeeeeeeeeeccecccccce 
c ‘angle.f' READS FILES 'fort.7' & ‘fort.3" WRITES OUTPUT TO FILE 'fort.10' 

c PROGRAM READS THE OUTPUT FROM RVCQ3D.f AND GRAPE.f AND GENERATES A 
c DATA FILE CONTAINING THE Y-POSITION AND INLET FLOW ANGLE AT A DESIRED 
c LOCATION UPSTREAM OF THE BLADE IN TERMS OF PERCENT CHORD. DATA IS 

c CALCULATED USING SIMPLE LINEAR INTERPOLATION. 

c SEPTEMBER 29, 1994 WRITTEN BY LT. DREW WILLIAMS 
eccecceccecececececeeeeceececeececeeeceeeeeceeeeeeeeeeeeeececeececeeeeececeeceeceeeceecccecece 


real beta(300) 

dimension x(300,300),y(300,300),q(300,300,4) 
dimension xx(300),yy(300),u(300,300),v(300,300) 
dimension uu(300),vv(300) 


open(10) 
rewind(10) 
pi=3.1415926 


c Read in the C-grid created by GRAPE 
read (7,*) ni,nj 
read (7,*) ((x(i,j),i=1 ni) j=1,n)), 
& ((yGij),= 1 ni) j=1,n)) 
read (7,*) mtl,mil 
¢  write(10,*) x(mtl,1) 


c Read in the Solution file created by RVCQ3D 
read (3,*) ni,nj 
read (3,*) fsmach,alpha,reyn,1ters 
read (3,*) (((q(iJj,K),1=1,mi), j=1,nj),k=1,4) 


c¢ Calculate the velocities u and v at every grid point 
do 100 1 = 1,ni 
do 200 j = 1,nj 
u(ij=atiJ,2)/qGij,1) 
v(i,j)=a(ij,3)/a(ij,1) 
200 continue 
100 continue 


c¢ Define the Chord length and then prompt for an upstream location 
chord=x(mtl, 1) 
write(*,*) ‘Input Upstream Chord Location in Percent (i.e., -10, -20,...)' 
c¢ (NOTE: You must be aware of the maximum upstream position available from your 


c GRAPE output. Beyond the maximum upstream location, the values calculated 
C are all the same as that at the maximum.) 

read(*,*) ax 

axperc=ax/100 

vertx=axperc*chord 
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Cc write(10,*) chord 
Cc write(10,*) vertx 


c Bracket the vertical x-position 
c (Note: Depending on your grid size, the ‘i’ index may have to be adjusted to 
Cc ensure the entire inlet flow region is covered.) 
do 300 1 = 1,240 
do 400 j = 1,nj 
flag=0 
if (x(i,j).gt.vertx) jsx1=] 
if (x(i,j).1t.vertx) then 
jsx2=j 
flag=2 
goto 500 
endif 
400 continue 
500 continue 
if (flag.gt.1) then 


¢ Calculate the u and v velocities and compute the inlet angle 
c Linear Interpolation in the form of: ynew = yold + (dy/dx)*dx 
yy(i)=y(i,jsx1)+(vertx-x(i,jsx1))*(y(i,jsx1)-y(i,jsx2))/ 
(x(i,jsx1)-x(i,jsx2)) 
uu(i)=q(1,jsx 1,2)+(vertx-x(i,jsx1))*(q(i,jsx1,2)-q(i,jsx2,2))/ 
& (x(i,jsx1)-x(i,jsx2)) 
Vv(i)=q(i,jSx1,3)+(vertx-x(1,jsx1))*(q(i,jsx1,3)-q(i,jsx2,3))/ 
& (x(i,jsx 1)-x(i,jsx2)) 
beta(i)=atan(vv(i)/uu(i))* 180/pi 
write(10,*) yy(1),beta(i) 
endif 
300 ~=— continue 
stop 
end 


100 


NG), 


HE . 
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